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Abstract 

We establish exact inequalities for the structure-function scaling exponents of a pas- 
sively advected scalar in both the inertial-convective and viscous-convective ranges. These 
inequalities involve the scaling exponents of the velocity structure functions and, in a refined 
form, an intermittency exponent of the convective-range scalar flux. They are valid for 3D 
Navier-Stokes turbulence and satisfied within errors by present experimental data. The in- 
equalities also hold for any "synthetic" turbulent velocity statistics with a finite correlation 
in time. We show that for time-correlation exponents of the velocity smaller than the "local 
turnover" exponent, the scalar spectral exponent is strictly less than that in Kraichnan's 
soluble "rapid-change" model with velocity delta-correlated in time. Our results include as 
a special case an exponent-inequality derived previously by Constantin & Procaccia [Nonlin- 
earity 7 1045 (1994)], but with a more direct proof. The inequalities in their simplest form 
follow from a Kolmogorov-type relation for the turbulent passive scalar valid in each space 
dimension d. Our improved inequalities are based upon a rigorous version of the refined 
similarity hypothesis for passive scalars. These are compared with the relations implied by 
"fusion rules" hypothesized for scalar gradients. 
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1 Introduction 



Much progress has been made recently in the understanding of anomalous scaling for the prob- 
lem of randomly advected scalars |Kr9|, pYGj [FGLP| , |GKl| , pGKj , |CFKL| , |CT| , |CFL| , ggl] , |SS2 |. 
The dynamical equation of the model is 



(dt + v(r, t). V r )0(r, t) = k A r 0(r, t) + /(r, *) 



(1) 



with 6(r,t) the scalar field, f(r,t) a (stochastic or deterministic) source, and v(r, t) a random 
incompressible velocity field. The issue of interest is the presumed scaling law 



S p (£) ~ 



(2) 



as 



-> for the scalar structure functions S p (£;6) = (\Ag9\ p ). For a high Reynolds number 
turbulence and for molecular diffusivity of the order of magnitude of the molecular viscosity, or 
greater, k > v, the dimensional theory of Obukhov flObu ] and Corrsin [Corr] implies that 



S p (£) ~ (x/e 1/3 ) p/2 Z p/3 



(3) 



in which \ = K ([^@] 2 ) is the (mean) scalar dissipation and e is the dissipation of kinetic energy. 
The scaling law Eq.(||) is supposed to hold for L ^> £ ^> (K 3 /e) 1//4 , where L is the length scale 
of the scalar source, assumed less than the integral scale of velocity. This specifies the inertial- 
convective range. Thus, Cp(^) = p/3 in the classical theory of this range. On the other hand, 
for k C v there is another range, L S> £ 3> (K 3 /e) 1 ^ 4 , in which it is now assumed that L is 
at, or smaller than, the Kolmogorov dissipation scale (z/ 3 /e) 1//4 . Over this range, the so-called 
viscous- convective range, the theory of Batchelor [ |Bat| ] proposes that 

(4) 



S p (e)~( x /< r y/*l g°*(£/L), 



with a a mean shear strength. Here, Cp(^) = formally for all p. Whereas [Obu, Corr| predicted 
Q p {0) to be a linear function of index p, it is now generally expected that these exponents are 
some nontrivial concave functions of p, i.e. that there is anomalous scaling [AHGA]. 



The new work on this problem cited above mostly deals with a special model in which the 
Eulerian velocity field is zero-mean Gaussian, delta-correlated in time 

(^(r,%(r',i')> = V i3 (v - r')5(t - t'). (5) 



This so-called "rapid-change model" was first investigated by Kraichnan [Kr68], who observed 



that for this case the infamous closure problem is absent: the iVth-order correlator of 9 obeys 
equations depending only upon itself and lower order correlators. The recent analytical investi- 
gations explore particular limiting regimes: the case with space-dimension d 3> 1 in JCFKL , [CI 



and the case with eddy-diffusivity exponent < ( <C 1 in |GK1 , BGKJ . The latter exponent is 



defined by the assumed scaling relation for the (Richardson) eddy-diffusivity tensor 
Kij{r) = \f ds {[v i (r,t)-v i {0,t))[v j {r,s)-v j (0,s)]) 

£ J -co 

= V ZJ (0)-V tJ (r), (6) 

that 

Ka (r) ~ D ■ (d - l)<%r C + D • Cr C - ^f) (7) 

for small r. Note that Eq.(^) is an analogue for this model of Taylor's 1921 exact formula for 
the eddy-diffusivity |Tay| (which involves instead Lagrangian velocities). 

It is our purpose here to consider the problem with a finite time-correlation of the convecting 
velocity field. This includes the realistic case where the velocities are turbulent solutions of 
Navier-Stokes dynamics. In addition, our results apply to a model recently considered | CFL|| 



with Eulerian velocity field taken as a Gaussian with covariance obeying dynamical scaling 
Vij(r,t) = ^{[v i (r,t)-v i (O,t)][v 3 (r,0)-v j (O,0)]) 



Dr< 



(8) 



with r r = Ti(r/L) z . Our results in the finite correlation-time models shall be applicable to 
both the limiting regimes studied for the rapid-change model: tf> 1 and < ( <C 1. However, 
rather than asymptotic formulae for the scaling exponents, we shall derive exact inequalities. 



One interest of our results is that they point up some significant differences between the zero 
and finite correlation-time problems. 

Our simplest set of inequalities are based upon the following relation: 

(A e v {l [A e 0] 2 ) = -^ X -i, (9) 

This equation is valid for L 3> £ ^> no, where L is the length scale of the scalar source and 
rjD is a dissipation length, given by the Obukhov-Corrsin length (K 3 /e) 1 / 4 |Obu, Corrf| for high 



Reynolds number Navier-Stokes turbulence and by (n/D) 1 ^ for the model of | CFI| (Eq.(||)) 



Eq.@ is a relation analogous to that of Kolmogorov for the third-order velocity structure 



function [K41-III| and it was proved for d = 3 by Yaglom in 1949 [Yag|. It simply expresses 
the constancy of scalar flux over the convective range of scales. By direct application of the 
Holder inequality, we will derive from this relation a basic set of inequalities relating the scaling 
exponents of pth structure functions of the scalar Cp($) with those of the velocity Cp( v )- The 
implications of these results will be discussed in Section 2. 

Furthermore, we shall derive an improved set of inequalities by means of a refined similarity 
relation (RSR) for the passive scalar. In a precise version stated below, the RSR we prove is 

^ amdia^dp (10) 



in which X^(r) is a local scalar flux to scales < t at space point r. (Cf. [ AHGAfl ). The 



corresponding inequalities will involve the intermittency exponent t p (X) of the scalar flux 

(ii) 

which measures the increasing spatial spottiness of the flux as I — > 0. These inequalities rigor- 
ously establish an intuitive fact: that convective-range intermittency of the scalar flux implies 
anomalous scaling of the scalar structure functions. These results are given in Section 3, along 
with some general discussion of refined similiarity hypotheses for passive scalars, including the 
relation to "fusion rules" proposed for scalar gradients. 
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2 Yaglom Relation Inequalities 

We shall sketch here very concisely the proof of the Yaglom relation for any space dimension d. 
See also [ jSKS [, and, for more details, | |UF| j and Appendix II of Eyl]. The first step is to define 



a (mean) "physical space scalar flux", via 

(12) 



X(l)=-~(9(r,t)9(r + l,t)) 



conv. ,i=0 

The subscript "conv." indicates that only the convective terms in Eq.(||) are used. A simple 
calculation using incompressibility and spatial homogeneity gives X(l) = — j Vi-(Aiv[Ai#] 2 }. 
Assuming also spatial isotropy, the vector A(l) = (Aiv[Ai#] 2 ) is (for d > 2) of the form 
A(l) = t4||(£)1, where 1 is the unit vector in the direction of 1. In the convective range of 
length-scales I with constant mean scalar flux, the Eq.(^) becomes —4% = Vi-A(l), or 



d - 1 . _ dA\\ 



±X = —^A ll (e) + -J-(e). (13) 



The only solution of this equation regular for £ —* is 

A ]l (£) = -^£. (14) 

This completes the derivation of the Eq.(^). It is useful to remark here that the Yaglom relation 
does not hold in the Kraichnan "rapid-change" model. In fact, the lefthand side of the Eq.(|]) is 
not even well-defined in Kraichnan's model, since it is expressed by a product of Agvu and [A^#] 2 
at a single instant. However, with the delta-correlation in time, the velocity is a distribution- 
valued (generalized) process and the single-time values are not defined. The analogous result 
for Kraichnan's model is Kr68(| 



ft W ~ zr^i^^ ( 15 ) 

which, like the Yaglom relation for the finite correlation-time velocity statistics, is an exact 
result in the delta-correlated model. 

A simple set of exponent-inequalities follow from the Yaglom relation as a straightforward 
application of the Holder inequality. The inequalities for C, p {9) involve as well the exponents 



Cp(y) of the (absolute) structure functions of velocity, S p (£;v) = (\A£v\ p ) ~ £(p( v ) . They are 
simplest to state in terms of the exponents a p {6) = ( p (0)/p, a p (v) = ( p (v)/p. By the Yaglom 
relation and the Holder inequality, 



4 X £ = (A e v {l [A e 9f 



(16) 



for 



2 1 
p>2, - + - = 1. 

P Q 



(17) 



As this relation must hold in the infinitely long convective range as i — > 0, it follows that 



2aJ9) + aJv) < 1. 



(18) 



We have used the isotropic form of the Yaglom relation, but this is inessential (see | Eylj| ). The 



special case of Eq.(18) for p = oo, q = 1 was previously derived by Constantin and Procaccia 
|CP1 | on the basis of estimates for Hausdorff dimensions of scalar level-sets: see Eq.(4.2) there. 
^ Our very simple derivation here shows that it belongs to a family of inequalities which are a 
consequence just of the constancy of scalar flux. These inequalities express a complementarity 
between the regularity of the velocity and scalar fields: if the (Besov) regularity exponent of 
velocity a q (v) is "big" then the corresponding scalar exponent a p (9) must be "small." 



lr The authors of [ |CP1 ] were not specific about which order p, q of exponents were involved, but inspection of 
the proof shows that p = oo and q — 1 — in our notation — was used. Actually, most of their proof generalizes to 
general p, q satisfying Eq.(^), except for the result on dimensions of level sets Eg, Dn(Eg) < d — <7oo(#), their 



Eq.(1.8), which used p = oo. Incidentally, we disagree with the conclusion of |CP1| that the estimate is sharp, 



i.e. an equality. This claim in Section 4 is based on an opposite inequality, 2a OB {6) + ai(v) > 1, supposed to 



be derived in [CP2], which we dispute. In fact, examination of the argument of [SP2| shows that an implicit 
assumption was made that a spatially uniform diffusive cutoff exists, below which length-scale the scalar graph is 
smooth. This assumption was used in the selection of ro in their inequality Eq.(4.23). To bound vold(G(B)) from 
below requires that G(B) is smooth on scales < ro so that a uniform choice of ro may be made. It is possible, 
however, that there is a "local fluctuating cutoff," as postulated in some multifractal pictures: see 
case, the uniform choice of ro would fail and the result, as well as its proof, might break down. 
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In that 
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The basic inequality Eq.(|l8|) may also be written as 



Cp(*)<f(l (19) 



for p > 2, - + - = 1. This relation may be considered under various special assumptions. If the 
scaling exponents of velocity v are taken to be K41, i.e. o~ q (v) = ^ for all q > 1, then it follows 
that 

C P (0) < I (20) 

for p > 2. Thus, the Obukhov-Corrsin predictions for the inertial-convective range appear as 
upper bounds. Likewise, if the velocity field is assumed smooth, or o~ q (v) = 1 for all q > 1, then 

(T P (6) < 0. (21) 

The smoothness assumption would hold, for example, for the velocity field in the viscous dis- 
sipation range, and then the Batchelor exponents for the viscous-convection range appear as 
upper bounds. More generally, if it is assumed that cr q (v) = h for all q > 2, then 

Cp(*)<f(l->0- (22) 

This assumption on the velocity scaling exponents corresponds to a "monofractal" field and 
would be true for convection by a Gaussian random velocity field. 

The last inequality has some interesting consequences for the model studied in [ |CFL ]. In 



that model the velocity field is spacetime Gaussian with covariance satisfying dynamical scaling, 
Eq.(|8|). Setting t = in that equation, it is easy to see that 2h = ( — z or 

C = 2h + z. (23) 

As pointed out above, £ has roughly the interpretation of an eddy-diffusivity exponent analogous 
to the Richardson exponent Cr = | [Rch ] . This is not entirely accurate since the exponent 



appears in the scaling law Eq.(||) hypothesized for Eulerian velocities in the model of [CFL|. 



For z < 1 Eq.(|8D is not an accurate representation of Eulerian time correlations, which will 



then be dominated by convective sweeping. Nevertheless, keeping to this terminology, there 
is also for fixed "eddy-diffusivity exponent" ( a complementarity between the magnitudes of 
velocity regularity exponent h and dynamical scaling exponent z: if one is big, the other is 
small. Of course, this is just due to the simple heuristic that eddy-diffusivity K% ~ v p t I- A 
similar relation should also hold for Navier-Stokes turbulence, i.e. £ = C2(v) + z, except that 
there the scaling exponent z will correspond to Lagrangian time-correlations. It is only for 
Lagrangian time functions that the dynamical scaling can hold and, furthermore, the exact 
Taylor formula Eq.(^) involves such correlations. By means of Eq.(^), the main inequality 
Eq.([l9|) may be reexpressed as 

C P (0)<~(2 + 7), (24) 
where 7 = 2 — £. In fact, Eq.(|2~3[) states that 



l-h = ^, (25) 

so that it is a direct consequence. Eq.(^) holds for all p > 2 in the model of | CFI| ; furthermore, 
it will hold also for p = 2 in Navier-Stokes turbulence, if C = Cz{ v ) + z as expected. 



In [CFL| an expansion about the "rapid-change model" was developed in the magnitude e 
of the correlation-time. They employed a particular choice of dynamical exponent 2 = 7. Their 
plausible physical argument for this choice ran as follows: the "local turnover time" of scalar 
eddies at scale £ is t% ~ £/v£ ~ £}~ h and this defines a "local turnover value" of the dynamical 
exponent z = 1 — h. Note by using Eq.(23) that this value is achieved precisely when ( = 1 + h, 
or 7 = 1 — h. In other words, the local turnover exponent is obtained when z = 7. For z > 7 
the velocity at vanishingly small scales changes randomly at a faster and faster rate relative 
to the evolution time of the scalar eddies. Hence, it is plausible that the predictions of the 
"rapid-change" model will hold in that case. On the other hand, for z < 7, the realizations 
of the velocity field are selected randomly at a rate which goes to zero compared to the scalar 
cascade rate, i.e. the velocity field randomness is "frozen in." We shall now observe that in the 
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latter case of quenched randomness, or z < 7, the scalar spectral exponent 

(2(0) < 7, (26) 
with strict inequality. This is direct from Eq. ( |24| ) . More generally, for p > 2 

Cp(0)<f-7- (27) 

Hence, all of the £ p 's for p > 2 are strictly smaller than the "classical" values. This is not so 
surprising for the higher p-values, p > 2, since this is a familiar situation usually associated 
to "intermittency." The present strict inequalities, including the unusual case p = 2, have 
a different origin. They arise just from the condition of constant mean flux, which requires 
smaller Cp($)' s when h is bigger. However, z smaller than the local turnover exponent at fixed 
value of "eddy-diffusivity exponent" £ requires h bigger than its classical value 1 — 7. Indeed, 
1— j < h < 1 — z from Eq.(^) if z < 7. 

Another important comparison between the "rapid-change" model and the finite time- 
correlation cases arises from the dimension-dependence of the Yaglom relation. It is clear 



from Eq.(15) for S2CO m the rapid-change model that it goes to a finite limit as d —* 00 if 
and only if D oc Dq/cP, with Dq fixed. In fact, it follows from the work of Chertkov et al. 



in [CFKL] that with that choice of (i-dependence of D all scalar correlations have a nontrivial 
limit as d — > 00, and, in fact, correspond to the correlations of a Gaussian field. Note that this 
dependence of D implies that each of the i = 1, ...,d components of the velocity field Vi(£) at 
scale I has typical magnitude and, likewise, each of the d 2 components of the strain tensor 
<T£ = i[(Vv^) + (Vv;) T ] are of order It may be argued on the basis of theory of random 



matrices that the typical strains along principal axes, the eigenvalues of 07, are order I 7 as 



d — > 00 : see [Kr74A, FPRj. This seems to be the correct scaling for a nontrivial limit, since the 



strain magnitudes gives the rate ~ 1/tg of scalar cascade. The Yaglom relation shows that 



2 Actually, the principal eigenvalue of <j£ may be expected to determine the rate of scalar cascade: see [Kx74A|. 
However, for random Wigner matrices with asymptotic semicircle distribution of eigenvalues, it is known also 
that the leading eigenvalue is within 0(l/d 2//3 ) of the right edge of the spectrum. 



the matter is not so simple for the finite time-correlation situation. For that relation a scaling 
v ~ 1/d is required to obtain a finite limit as d — > oo. This does not contradict the results 
of |CFLj| , since they take e = Tg/tg as the small parameter of their expansion. Working this 
through, one finds that this amounts to taking vg ~ -^L^ 1- ^ and Tg ~ e • For any d, this 
correctly recovers the delta-correlated model in the limit e —* 0. In fact, e is just the quantity 



denoted t* in [Kr6q|. For the validity of the Yaglom relation it is therefore required that e ~ d 



and that is clearly incompatible with the condition of [CFL] that e <C 1 at large d. [] 



3 Refined Similarity Inequalities 

We shall now derive inequalities which improve those from the Yaglom relation. The basic idea 
of the proof is a scaling relation between the local scalar flux variable Xg(r) and the difference 
variables of velocity and scalar at the same point r, Eq.(|l0|). This is an analogue of the refined 



similarity hypothesis (RSH) in 3D which — in the version of Kraichnan [Kr74B] — states that 
local energy flux scales as n^r) ~ [Agv(r)] 3 /£ in terms of the velocity difference at the same 
point. The proofs given below follow closely methods used in our discussion of the 3D RSH in 



[Ey2| and the 2D RSH for vorticity scaling exponents in [Eyl]. If we assume 



(\X e \P) (28) 



then there follow heuristically from Eq.(|10|) relations between the exponents Cp(0)>Cg( w ) an d 

T r (X). The exponents t p (X) measure the increasing spatial intermittency or "spottiness" of the 

scalar flux at decreasing length-scales. In fact, since (Xg) = x over the long convective interval 

of £, it may be expected that the (concave in p) exponent t p (X) is negative for monents p > 1. 

The corresponding growth in moments of Xg as t — > reflects the increase in its fluctuations. As 
3 The proportionality of single-time velocity realizations to e -1 / 2 here, as well as in Kraichnan's original 1968 
derivation, makes clear that such single-time values do not exist in the delta-correlated model obtained by the 
limit e — > 0. Only integrals over some finite time-interval are well-defined, and the lefthand side of the Yaglom 
relation, for the idealized white-noise limit, is a meaningless expression. 
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we shall establish below, the intermittency of the convective range scalar flux implies anomalous 
scaling of the scalar structure functions over that same interval. It is this connection between 
intermittency and anomalous scaling which is the essence of Kraichnan's RSH [Kr74B]. After 



these results are derived as theorems below, we shall comment on the relationship with other 
refined similarity hypotheses for passive scalars recently proposed [Hos, ZAH, |5KS |. These 



latter hypotheses are motivated by the original Kolmogorov RSH |K62 | which involves space- 
averaged dissipation rather than local flux. 

We first must introduce an appropriate definition of the local scalar flux. It is most easily 
done using a smooth filtering technique to differentiate the large-scale and small-scale modes. 
This is the same method used in the large-eddy simulation (LES) modelling scheme and in our 
earlier discussion of the 3D case [Ey2|. Here we apply the filter to the scalar equation Eq.(|l]). 



That is, we consider the "large-scale scalar field" defined as the convolution field 9e = Gi * 9, 
with some suitable filter function Gg. The resulting equation is 

d t 9 £ (r,t) + V-(v e (r,t)e e (r,t) +] e (r,t)) = 0. (29) 



The large-scale velocity field is likewise defined by = Gi * v. Note that jg = (v9) e — V£#£ is a 
space-flux of the scalar induced by the turbulent convection (eddy diffusion) . A main ingredient 
of our proofs is the following exact formula for this turbulent flux: 

j £ (r, t) = [A0(r, t)Av(r, t)] t - [A0(r, %[Av(r, t)] e . (30) 

Here [f]e = J d 2 s G^(s)/(s) is the average over the separation-vector s in the difference-operator 
A s with respect to the filter function Ge(s). See [ |Eyl , Ey2|. This relation shows that j^(r) ~ 
A^v(r) • A^#(r) at each space-point r. 

Recall that the scalar-intensity integral K(t) = \ Ja ^ 2 (^) ^ s formally conserved by the full 
dynamics. From the Eq.(]29|) for the large-scale scalar field it is straightforward to derive by the 
standard methods of nonequilibrium thermodynamics a local balance equation for its large-scale 
intensity Kg = It has the form 

D t K e (r, t) + V-D^r, t) = -X e (r, t). (31) 
11 



Here, D t represents dt + v^(r, f)-V r , the convective derivative associated to the large-scale 
velocity, 

D*(r,t) =e e (r,t)j e (r,t) (32) 

represents space-transport of the large-scale intensity by convective diffusion, and the scalar 
flux 

X e (r,t) = -Ve e (r,t)-j e (r,t) (33) 
represents the scalar transfer to the small-scale modes. In a homogeneous, stationary ensemble 



the lefthand side of Eq.(31) has vanishing average. In a steady-state with constant mean flux 



X of scalar substance to high wavenumbers, the average {Xg) = x, a constant, for I lying in 



the convective interval. Together with Eq. ([}(]), the formula Eq.(32) for scalar flux shows that 
X e (r) ~ A £ v(r) • [A £ 0(r)] 2 /£, which is the RSR, Eq.flHJ). It is the exact equations Eqs.@),(||) 



which are the precise form of our RSR, applicable even without assumptions of local isotropy 
or other statistical properties. They are essentially kinematic in nature, based only upon the 
conservation properties of the underlying dynamics. 

From these exact relations, there follow refinements of the previous exponent-inequalities. 
In fact, it follows from the (generalized) Holder inequality that for 



the ordering holds that 



2 11 

- + - = -, (34) 
pgr 



fr(X)/r ^ (\x e \ r ) 1/r 

{\A t 0\P) 2 /P{\A e v\*)Vi 



< 



and thus, 



2a p (9) + a q (v) < 1 + Iri^l, (36) 



This is our main result on the scaling exponents. The derivation requires only the exact 
kinematical relations, Eqs.(|30|),(|33|), rather than the heuristic form of the RSR, Eq.(p!o|). For 
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the details of this, see the Appendix of [Ey2]. Note that r = 1 recovers the previous inequalities, 
Eq.(|l§|), if t\(X) = 0. However, as noted above, intermittency of scalar flux will imply that 
T r (X) < for r > 1, and then the inequalities are sharpened. This is particularly easy to see 
for the case of a "synthetic" turbulent convection by a Gaussian random velocity field. For 
the Gaussian field (or any "monofractal" field) cr q (v) = h for all q. Thus, taking q — > oo in the 
above inequality, one easily obtains r = p/2 and 



C P (e)< F -{i-h) + T P _{x). 



(37) 



Thus, Cp($) will be strictly less than the "classical exponent" C; lass = |(1— h) and the "anomaly" 
is exactly an intermittency exponent of the convective-range flux. This result implies that 
convection by a regular-scaling random velocity will nonetheless lead to anomalous scaling for 
the scalar it passively convects, if the associated flux variable develops strong fluctuations. 
It is worthwhile to make a comparison of these results with the other recently-proposed 



RSH for passive scalars [Hos, ZAH, SKS]. While our formulation is motivated by the 1974 
"revisionist" RSH of Kraichnan, involving flux, the RSH explored by the above authors is an 
adaptation of that originally proposed by Kolmogorov | K62f| . That is, it is hypothesized that 
the random variables 



V&(r,l) = A,0(r; 



(&<to) 1/6 



(38) 



(^(r)) 1 / 2 ' 

defined in terms of volume-averaged dissipations £i(r) and Xe( r ) °f velocity and scalar intensi- 
ties, respectively, have conditional distributions given values £g and Xii which are independent 
of the local Reynolds number Reg and local Peclet number Peg when those are both large. 
That is, the variable Vq is supposed to have a universal distribution in the inertial-convective 
range of t. If this relation is combined with the original Kolmogorov RSH, then it is easy to 
infer likewise the existence of a random variable Wg, universal in the same sense, such that 

Atv ■ [A e 9] 2 - W e ■ Xlt (39) 

See [ [SKS , |Z AH ] . Given the Kolmogorov RSH for velocity-differences, this last relation is, in 
fact, equivalent to the RSH for passive scalars proposed in Hos, ZAH, 3KSJ. It provides a 
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natural bridge between that RSH for passive scalars and the one established here. Since we 
have shown that X( ~ A?v ■ [A£0] 2 /£, the above relation may be more or less paraphrased as 

Xt~Wg- XI- (40) 

In other words, the ratio Xi/xi — Wg is a random variable whose distributions conditioned on 
fixed ££, xt are universal in the inertial-convective interval of I. Again, this is essentially just a 



reformulation of the RSH of [Hos, ZAH]. If it holds, then a simple consequence is that 



(\X e \P) = {\W e \ p \e t ,xe,e)-{x P e ), (41) 

and the coefficient (|I^g| p |£^, Xii^) ls -> 111 the inertial-convective range, just a constant factor 
Therefore, in particular, t p (X) = t p (x) for all p, and the intermittency-exponents of 
the convective-range flux X^ and the scalar dissipation, volume-averaged over the same length- 
scales, Xlt are the same. In that case, all of the inequalities previously rigorously derived in 
terms of t p (X) hold also for T p (x)- In terms of providing a theoretical foundation to the RSH 
for passive scalars, it may be easier to proceed by starting with Eq.(flc|). 

Another interesting comparison involves the additive fusion rule (AFR), which was proposed 



originally for the turbulent velocity-gradients | Ey3 , LLJ| . Recently, the straightforward extension 



of these rules to the scalar-gradients has received some analytical support in Kraichnan's "rapid- 



change" model [FGLP, BGK, CFKL, CF|. As we now explain, it happens that in this model 
the AFR and the RSH lead to identical relations between scaling exponents. The AFR states, 
in a schematic form, that 

[(V6) P1 ] ■ ■ • [(V0) p «] ~ [(V0) Pl+ - +p "]. (42) 

We ignore for the sake of this argument the vector character of the scalar-gradient which, prop- 
erly, should be taken into account (Cf. [ Ey3| ] ) . The quantities [(W) p ] are so-called "renormal- 



ized composite variables." This means simply that one defines them as the limit of pth-powers of 
scalar-gradients at the same space-point, in the model with r]£> > 0, but multiplicatively renor- 
malized by an appropriate power of r]r>- After nondimensionalizing the variable (according to 
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its "canonical" or engineering dimension) there may still be required a factor Z(rj£>) ~ (^lO p 
to make the correlations finite in the limit tjd 0. The exponent x p is the so-called anomalous 
scaling dimension of the variable [(V#) p ]. Observe, in this context, that the scalar dissipation 
x(r) = k|V#(t)| 2 — when divided by its mean \ — is nothing more than [(V#) 2 ]. The precise 
meaning of the schematic result Eq.fli^) is that, inserted in arbitrary correlators at separated 
points, 

[(V0r](A- ri )...[(V^n(A-r n ) 

^ \x P1 +-+Pn-x Pl *pn[(V60 Pl+ "' +p "](O), (43) 



in the limit as A — > 0. It is easy to show, as in [Ey3], that the above AFR leads to a "multiscaling 



law" for scalar structure functions (now without absolute values ) of the form 

([A e 9] p ) ~ e0-- x i)p+ x p . (44) 

Moreover, if the short-distance expansion is applied to the moments of the volume-averaged 
dissipation, it is immediately obtained that 

(x P e ) ~£*»p-p*'. (45) 

In Kraichnan's model, x\ = 1 — (7/2) and X2 = 0. Thus, Eq.(p3|) leads to ( p (0) = p(j/2) + x p 
and Eq.(^) yields T p (x) = x 2p- Hence, they together give 

( p (0)=P-~ + tz(x). (46) 



This should be compared with the result of RSH for the model, Eq.(37) [taken as an equality 
and with r p (x) replacing r p (X).] Clearly, they are the same. It was already shown in [|CF| that 
the RSH holds in the white-noise model. Our point here is that the RSH is a consequence just 
of the AFR. 

4 This happy situation does not, however, hold for the original AFR applied to velocity-gradients, nor even 
necessarily for scalar-gradients in true turbulence. In the case of velocity-gradients, the application of the AFR 
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One use of our exact inequalities is as a check on experimental data for scaling expo- 
nents, since these employ a variety of assumptions and approximations (Taylor hypothesis, 
one-dimensional surrogates, etc.) For that purpose, we may cite the experimental results of 
[ AHGA|1 

C 2 (0)«O.65, C3(0)~O.82, C4(0)~O.95 (49) 



for the scalar exponents, and [HvW] 



a 2 (v) w 0.355, a 3 (v) = 0.333, a^v) < a 17 (y) rj 0.211 (50) 

for the velocity exponents. We shall here assume crn(v) « cr oc(v) since the graph of ( p (v) is 
close to linear for p of order 20. In that case, if we make comparison, for simplicity, with the 
Yaglom inequalities, we find that 

0.65 - ( 2 (0) < 1 - a^v) ^ 0.78 (51) 

O.82«Cs(0)< -(l-<r 3 (v))« 1.00 (52) 

0.95 « ( 4 (6) < 2(1 - a 2 (v)) « 1.29 (53) 

All of these inequalities are well-satisfied by the data. The fact that there is a considerable mar- 
gin between the upper and lower limits is also consistent with an intermittency correction from 
the scalar flux. Determination of the latter from DNS or, experimentally, from the surrogate 
Agv [A£0] 2 /£, would be of interest. 



analogous to the above leads to the relation Ey3] 



whereas the RSH leads instead to 



&(«)=P'^+r 5 (e), (47) 



Cp(«) = |+r f (e). (48) 



These are not equivalent and the RSH result seems to be in better agreement with the experimental data 
(Sreenivasan, private communication, 1995). The failure of the AFR in this context could be related to the 



existence of a hierarchy of viscous cutoffs for the velocity field FV , Ey3 , or, possibly, to a naive disregard of the 
tensorial character of products of velocity-gradients. 
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